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Abstract 

We study both theoretically and experimentally anchoring properties of photoaligning azo-dye 
films in contact with a nematic liquid crystal depending on photo-induced ordering of azo-dye 
molecules. In the mean field approximation, we found that the bare surface anchoring energy 
linearly depends on the azo-dye order parameter and the azimuthal anchoring strength decays 
to zero in the limit of vanishing photo-induced ordering. From the absorption dichroism spectra 
measured in the azo-dye films that are prepared from the azo-dye derivative with polymerizable 
terminal groups (SDA-2) we obtain dependence of the dichroic ratio on the irradiation dose. We also 
measure the polar and azimuthal anchoring strengths in nematic liquid crystal (NLC) cells aligned 
by the azo-dye films and derive the anchoring strengths as functions of the dichroic ratio which is 
proportional to the photo-induced order parameter. Though linear fitting of the experimental data 
for both anchoring strengths gives reasonably well results, it, in contradiction with the theory, 
predicts vanishing of the azimuthal anchoring strength at certain non-zero value of the azo-dye 
order parameter. By using a simple phenomenological model we show that this discrepancy can 
be attributed to the difference between the surface and bulk order parameters in the films. The 
measured polar anchoring energy is found to be an order of magnitude higher than the azimuthal 
strength. Our theory suggests that the quadrupole term of the spherical harmonics expansion for 
the azo-dye - NLC intermolecular potential might be of importance for the understanding of this 
difference. 
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I. INTRODUCTION 



When a nematic liquid crystal (NLC) is brought into contact with an anisotropic sub- 
strate, the energy of the NLC molecules in the interfacial layer and thus the surface tension 
[the excess free energy per unit area] will be orient at ionally dependent. The anisotropic part 
of the surface tension — the so-called anchoring energy — gives rise to the phenomenon 
known as anchoring, i.e., surface induced alignment of the nematic director along the vector 
of preferential orientation referred to as the easy axis. 

Over the past few decades anchoring properties of NLCs have been the subject of intense 
studies for both technological and more fundamental reasons. There are a number of surface 
ordering and anchoring transitions that were observed experimentally and were studied using 
different theoretical approaches [see, e.g., Refs. [110,0] for reviews]. 

Technologically, producing substrates with anisotropic anchoring properties is of vital 
importance in the fabrication of liquid crystal electrooptic devices. The traditional technique 
widely used to align liquid crystal display cells involves mechanical rubbing of ahgning layers. 
This method, however, has the well known difficulties related to physical damage, impurities, 
dust contamination and generation of electrostatic charge |j^]. 

An alternative photoalignment technique avoiding the drawbacks of the mechanical sur- 
face treatment was suggested in Refs. 0111,0. It uses linearly polarized ultraviolet (UV) 
light to induce anisotropy of the angular distribution of molecules in a photosensitive film . 

The photoalignment has been extensively studied in a number of different polymer svs- 
tems such as dye doped polymer laye rs t5 >, ^] , ci nnamate polymer derivatives ^, 2h 
and side chain azopolymers [3, Ej, llfit I1(tI. IitI . Eii Eoj ] . Light induced ordering in the pho- 
tosensitive materials, though not being understood very well, can occur by a variety of 
photochemically induced processes. These typically may involve such transformations as 
photoisomerization, crosslinking, photodimerization and photodecomposition (a recent re- 
view can be found in Ref. [20]). 

In this paper we examine anchoring properties of the films containing photochemically 
stable azo dye structures that were recently studied as new photoaligning materials for NLC 
cells (2l[l2^. Dependence of the surface anchoring strengths on the photoinduced anisotropy 
will be of our primary interest. 

More specifically, we are aimed to study the effects of the photoinduced ordering in 
azo-dye films on the polar and azimuthal anchoring energies. The key point is that the 
photoalignment technique provides a means for controlling the photoinduced ordering that 
affects anchoring properties of photoaligning layers by changing ordering of azo-dye molecules 
at the surface and, thus, the surface anchoring strengths. 

Recently, the anchoring properties of aligning photopolymer layers in relation to the 
photoinduced ordering were studied experimentally in Ref. j23|. The relationship between 
the rubbing strength and the azimuthal anchoring energy was discussed in Ref. 124] 

The photopolymer-NLC interface was also described theoretically in Refs. [2^ |2a] using 
a modified version of the variational mean field approach which is also known as the Maier- 
Saupe theory. By contrast, the azo-dye films have not yet received a proper attention and 
we intend to fill in the gap. 

The pap er is organized as follows. In Sec. |n] we apply the mean field theoretical ap- 
proach H Q to express the surface anchoring energy in terms of the tensorial order 
parameters which characterize angular distribution of the azo-dye and NLC molecules at 
the interfacial boundary surface. The general result is then used to derive the expressions 
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for the azimuthal and polar anchoring strengths that, in addition to the order parameters, 
depend on the harmonics of the intermolecular potentials. 

Experimental details are given in Sec. IIIIl The polymerizable azo-dye monomer SDA-2 
was used to prepare the photoaligning layers. Absorption dichroism spectra were measured in 
the films irradiated with linearly polarized UV light at various irradiation doses. Anchoring 
energy measurements were performed in NLC cells where NLC is sandwiched between the 
glass plates coated with the azo-dye film. 

In Sec. II VI we present the experimental results and apply the theory of Sec. |n]to interpret 
the data. Discussion and concluding remarks are given in Sec. Details on some technical 
results are relegated to appendices IXllBl 



II. THEORY 

In this section we begin with introducing general notations and apply the mean-field 
approach to express the Landau-de Gennes surface free energy in terms of both azo-dye and 
NLC order parameters. Expressions for the azimuthal and polar anchoring strengths, 
and We, are then derived from the orientationally dependent part of the surface energy in 
Sec. Ill (^1 In the concluding part of this section we consider effects of spatial variations of 
the azo-dye order parameter using a simple model formulated in Sec. IIIDI 



A. Order parameter and dichroic ratio 

Assuming that the unit vector, u = (sin ^ cos 0, sin 6' sin 0, cos 6'), directed along the long 
molecular axis defines orientation of a molecule in both azo-dye film and NLC cell, quadrupo- 
lar orientational ordering of the molecules can be characterized using the traceless symmetric 
second-rank tensor ||30|] 

Q(u) = (3u® u-I)/2, (1) 

where I is the identity matrix. The dyadic (0) averaged over orientation of molecules with 
the one-particle distribution function pa{r,u), describing the orient at ion- density profile of 
azo-dye {a = A) and NLC {a = N) molecules, is proportional to the order parameter tensor 



Po(r, u)Q(u)du = Po(r)S„(r), (2) 

where dn = sin6'd6'd0, pa(r, u) = pa(r)/a(r, u), Pa(r) = f pa(r, u)du is the density profile 
and /Q,(r, u) is the normalized angular distribution. The general expression for the order 
parameter is given in appendix^ [see Eq. ()A15|) ] along with technical details on the technique 
of irreducible tensors. 

Now we dwell briefly on the relation between the order parameter S\ characterizing 
orientational distribution of azo-dye molecules /a(u) and the absorption dichroic ratio 

where D\\ [D^] is the absorption coefficient measured for a testing beam linearly polarized 
parallel [perpendicular] to the polarization vector of the activating UV light which is directed 
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along the x axis, Eex = -Eex x. We shall also assume that the testing and the pumping waves 
are both propagating along the z axis which is normal to the film substrate. 

When the absorption tensor of an azo-dye molecule is uniaxially anisotropic with aij (u) = 
o'±Sij + (cr|| — a±)uiUj, its orientational average takes the following matrix form 

{a) = (aavI + 2AorSA)/3, (4) 

CTav = 0-|| + 2o-_L, Aa = 0-|| - (T^, (5) 

where the angular brackets (...) denote orientational averaging. 

In the low concentration approximation, the optical densities D\\ and D± are proportional 
to the corresponding components of the tensor ^ 

^11 ocpAKv + 2Aa5(t))/3, (6) 
I^±ocpAKv + 2Aa4^))/3, (7) 

so that the average absorption coefficient Dav is given by 

= + 2D^ oc PA (^av + 2/3 Aa [^^^ - ] ) . (8) 

When the absorption coefficient does not depend on irradiation dose (and, thus, on the 
order parameter), from the expression (jH)) we may conclude that anisotropy of the azo-dye 
film is uniaxial and Syy^ = Si^^ = = — 5'a/2. In this case we have 

Sa = 5a(3x®x-I)/2, R=—Sa. (9) 

0"av 

As is seen from Eq. (jH}, the dichroic ratio equals the order parameter only in the limiting 
case where absorption of waves propagating along the long molecular axis is negligibly small 
and o"_L 0. 



B. Anisotropic part of surface energy in the mean-field approximation 

In the previous section it was shown that the light induced ordering of azo-dye molecules 
can be described by the order parameter which is expected to affect the surface free energy 
at the nematic-substrate interface. So, in this section, the order parameter dependent part 
of the surface energy will be of our primary concern. 

In the case of a flat structureless substrate, the expression for the surface energy was 
originally obtained by Sen and Sullivan in Ref. l27l . Subsequently, similar results have been 



derived by using the mean-fleld approximation j28[ and the density functional theory [31 



Similarly to Ref. j2^] , we adopt the mean-fleld approach and use the Fowler approxima- 
tion for the one-particle distribution functions 

pt^{r,u) = H{z)pti{z,u), pA{r,u) = H{-z)pa{z,u), (10) 

where H{z) is the Heaviside step function which equals unity when z is positive and vanishes 
otherwise. 

Applying the mean-fleld theory 0] gives the Landau-de Gennes surface free energy as 
an excess Helmholtz free energy per unit area that depends on two pair intermolecular 
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potentials: (a) the potential of interaction between NLC molecules, f/N-N(i"i2, ui, U2); and 
(b) the potential of interaction between NLC and azo-dye molecules, f/A-N(i"i2) ui, U2), where 
ri2 = I"! ~ 1*2 is the vector of intermolecular separation and Uj is the orientation coordinates 
of the interacting molecules. The resulting expression is given by 

/O /'oo f- 

dzi j dz2 / duidu2 
■00 Jo J 



X 



Pa{zi, Ui)Va{zi2, Ul, U2)Pn(^2, U2) 
-^Pt^{Zi,Ui)V^{Zi2,Ul,U2)pt^{z2,U2) , (H) 

K.(^i2,Ui,U2) = / f/a_N(ri2,Ui,U2)da;i2d?/i2, (12) 

J A 

where A is the area of the substrate and Va is the potential averaged over in-plane coordi- 
nates. 

It should be noted that the potentials ?7n-n and J/a-n actually represent the perturbative 
part of interaction that can be treated in the mean-field approximation. They can be written 
in the form of expansion over spherical harmonics given in Eq. ()B1|) of appendix El For our 
purposes, however, it is more convenient to use the tensorial representation for the averaged 
potentials that was introduced in Ref. js^]. In appendix El the coefficients that enter 
this representation [see Eq. ()B4|1 ] are related to the coefficients, Vjj^j{z) with ji < 4, in the 



spherical harmonics expansion (jB4j). This relation is given by Eqs. ()B9jl - ()B12|l . 

Substituting the representation ()B4j) into Eq. (fTTjl and assuming homogeneity of the 
contacting phases, we obtain the Landau-de Gennes expression for the surface free energy 
in the final form: 

/5(Sn, Sa) = /n(Sn) + +/a(Sn, Sa), (13) 

/n(Sn) = Co z ■ Sn ■ z + Cn-* Tr(SN) 

+ cg)z.S^z + cgM2-SN-z]% (14) 

/a(Sn, Sa) = cJJ^ Tr(SN Sa) 

+ c^^ z ■ Sn Sa ■ z + cf' [z • Sn • z] • [z ■ Sa ■ z] , (15) 



where the coefficients are given by 



co = bf^-b^^\ (16) 

c^ = ^^ c« = -6«/2, (17) 

POO 

b^^=PaPN zPji\z)dz, (18) 



I3a\z) denote the coefficients in the representation ()B4|1 for the potential (fT^ . 

Eqs. f|T!?|l - (fTH|) can be viewed as a generalization of the expression by Sen and Sullivan |2j 
supplemented with the term /a(Sn, Sa) resulting from the interaction between NLC and 
azo-dye molecules. Note that this result can also be derived by constructing invariants from 
the order parameter tensors S^ and the normal to the substrate z. In this case, the surface 
of the azo-dye aligning film is treated phenomenologically as a bounding surface which, in 
addition to the normal z, is characterized by the order parameter Sa. 
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C. Bare anchoring energy 

Separating out the director dependent part of the surface free energy requires the order 
parameters of azo-dye and NLC molecules be substituted into Eqs. (|T3|l - (|T5|l . Since the order 
parameter at the surface may differ from its value in the bulk, we generalize the expression 
for the azo-dye order parameter as follows 

2Sa|.=o = 5a(3x ® X - I) + Pa(z ® z - y ® y). (19) 

Similarly, for NLC order parameter tensor at 2; = 0, from Eq. ()A15j) we have 

2Sn|^=o = 5'(3n (g) n - I) + P(m (g) m - 1 (g) 1) 

= (35 + P)fi® n + 2Pm® m- (5 + P)I, (20) 

where n = (sin O cos $, sin G sin $, cos G) is the NLC director, n ± m = cos7e^.(n) — 
sin7ey(n), e^^h) = (cos G cos $, cos G sin $, — sin G), ey{h) = (— sin $, cos $, 0) and 1 = 
n X m. 

Eqs. (fT^ and (pUj) suggest that the order parameters of azo-dye and NLC molecules 
though being both uniaxial in the bulk can be biaxial at the surface. In addition, the scalar 
order parameters Sa and S at the surface may also deviate from their values in the bulk. 

The surface free energy can now be expressed as a sum of two contributions 

/5(S,Sa) =W^(n,m)+/scai, (21) 

where W{h, m) is the orient at ionally dependent part of the surface energy. This part can be 
calculated by substituting the order parameters (|T9|l and pO|l into the surface energy (|T3| 
to yield the expression for the bare anchoring energy 

W{h, rii) =iV^(n-z)^ + M^(m-z)^ 
+ Ar,,(n.x)' + M,(m-x)' 

+ 4'V4 [(35 + P) (fi ■ z) ' + 2P (m • z) '1 ' (22) 
with the coefficients defined by the relations: 

4iV, = {3S + P)[q + {S - P)] , (23) 

4N, = c'i\3SA + Pa) (35 + P), (24) 

2M, = P[q - 24'^ 5], 2M, = c^I\3Sa + Pa)P, (25) 

q = 2co- {c^l^ + ) {Sa - Pa) + 2ci')pA - 2c^^\s + P). (26) 

The second term on the right hand side of Eq. ()2H) 

4/scai = - [2co - {Sa - PA)(cf + cf ) + 3c^I\Sa + Pa)] 
x{S + P) + (cSJ) + 4^){S + Pf + 24\3S' + P2). (27) 

is a quadratic function of NLC scalar order parameter S and the biaxiality P. From Eq. ()27p 
it, similarly to the anchoring energy ()22|) . depends linearly on the azo-dye parameters Sa 
and Pa. 
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For the anchoring energy ()22j). we consider the simplest case which occurs when the 
surface induced NLC biaxiahty P is neghgibly small and the quadrupolar term ^224 in the 
expansion of the intermolecular potential Vn can be ignored. Under these circumstances, 
setting P = = and = = 0, we arrive at the simplified formula for the anchoring 
energy 

iy(n) = iV,(n.z)' + iV,(fi.x)' (28) 
which agrees with the expression for the anchoring energy recently proposed in Refs. [ssi 

MM- 

From Eq. (j28|) it is clear that the easy axis is directed along the y axis, = y, only if the 



coefficients and Nz are both positive. In this case the polar and the azimuthal anchoring 
strengths. We and W^, are given by 

We = 2iV„ W^ = (29) 

From Eq. we immediately deduce a more explicit expression for the azimuthal an- 
choring strength 

W^ = w^[S^^^U=o-Si^^U=o], (30) 
2w^ = -3c^l^SU=o, (31) 

where notations indicate the plane z = as a surface separating the phases. 
Similarly, Eq. ()23|) gives the polar anchoring strength in the explicit form 

yy 

,(1) _ , (3h 



We = + ^^'^^if l.=o - [5(^)^0 - 5^^) Uo] , (32) 
2t.« = 3(c(J)+c(JVl-o. (33) 



The formulas (pn|l - (jHHjl will be subsequently used in Sec. HYIto interpret the experimental 
data. At this stage, it is worth noting that, for the order parameter Q, the relations 
and (jHHj) provide the inequalities 

ci')<0, cf+ci')>0 (34) 

as conditions for the anchoring strengths to increase linearly as the scalar order parameter 
Sa decreases. 



D. Model of spatially varying order parameter 

As it was pointed out at the beginning of the previous section, the surface order parameter 
tensor of an azo-dye film may differ from the bulk order parameter of the film . The 
latter, according to the experimental results presented in the subsequent section IIIH is 
uniaxially anisotropic with the in-plane anisotropy axis that is normal to the polarization 
vector of the activating UV light. In addition, the light induced scalar order parameter, 
which is proportional to the dichroic ratio 0, turns out to be negative, 5*^^'' < 0. 

From the other hand, assuming that the anisotropic part of the surface energy can be 
taken in the general form by Sen and Sullivan |2^, the boundary conditions may favor either 
homeotropic or planar alignment of the azo-dye molecules, thus, counteracting the action of 
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FIG. 1: Critical order parameter Sc as a function of the coupling constant g\ at various values of 
the coefficient g2- 



light. So, it can be expected that the effects caused by interplay between the light induced 
and the surface ordering are of importance in explaining the order parameter dependencies 
of the polar and azimuthal anchoring energies. 

In this section we discuss these effects on the basis of a simple phenomenological model 
formulated by using the polar representation ()A21|) for the azo-dye order parameter. The 
latter can be conveniently rewritten in the form 



= — sa cos^/), Pa = — V^SAsin?/', 



(35) 



where the angle tl) is shifted by vr so as to have the angle tl) vanishing in the bulk. 

In what follows we shall assume that, similarly to nematic liquid crystals [3 
the amplitude sa varies in space much slower than the angle tl). So, in our model, the 
amplitude will be fixed at its bulk value, sa = and we consider the limiting case of 

thick films in which the characteristic length of spatial variations of the angle il) is much 
shorter than the film thickness. In this case the film can be regarded as a semi-infinite 
sample filling the upper half space, z > 0. 

Technically, our task will be to find the spatially varying angle as a function of z that 
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FIG. 2: Components of the order parameter at the surface as functions of the bulk order parameter 
SA- Two cases are shown: (a) gi = 0.05 and g2 = 0.0; (b) gi = —0.2 and g2 = 0.4. 



minimizes the excess free energy per unit area, Fa, taken in the following nematic-like form 



Lsl{d,^Y + Bsl{l - cos(3^)) 



dz 



+ Gi sa cos('?/'o + vr/3) + G2 s\ cos^(-?/^o + 7i"/3), 



(36) 



where ipo = 4'\z=o and sacos('?/'o + 7r/3) = z ■ Sa ■ z\z=o. The first part of the excess free 
energy ()3fi|) is of integral form with the integrand describing the energy costs for deviations 
of ip from the equilibrium value, ip = 0. The gradient term of the energy density is taken 
to be proportional to ((9^Sa)^, whereas the other term gives an increase in energy caused 
by spatially uniform changes in ip. This term is written as a linear function of the angle 
dependent invariant ()A19|) . 4Tr[S^] = 3s|^ cos(3'?/^). 

For the order parameter (|T9|) . the surface part of the energy ()36|1 can be represented by a 
quadratic polynomial of z ■ Sa ■ z. By contrast to the elastic constant L and the coefficient 
B, the surface coupling constants, Gi and G2, can generally be negative leading to different 
boundary conditions. For example, if G2 = 0, minimizing the surface term requires the 
z-component of the order parameter, z ■ Sa ■ z = Si'^^ , to attain its maximal (minimal) value 
at the surface provided the coefficient Gi is negative (positive). These can be referred to as 
the homeotropic (planar) boundary conditions. 
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FIG. 3: (a) Scalar order parameter and (b) biaxiality in relation to the ratio z/^. 



The Euler-Lagrange equation for the free energy functional f)36p can be readily solved to 
yield the relation 



tan(3V^/4) = tan(3?/'o/4) exp(-z/0, 



(37) 



where 9^^ = 2L{Bsa) ^- This relation can now be substituted into Eq. (jHUj) to derive the 
free energy as a function of the angle ipQ. The result is 



.5/2 



sin^(3V'o/4) + giSA cos(V^o + 7r/3) 



+ g2si cos^(V'o + 7r/3), 



(38) 



where h = 4(2i?L)^/^/3 and Qi = Gi/h. The angle ^jj at the surface then can be found as the 
value of i/jq that minimizes the function (jHH|l on the interval ranged from —2tt/3 to 27r/3. 

Qualitatively, dependence of ipo on the coupling constant gi can be analyzed using ele- 
mentary methods. For g2 > 0, we find that the angle ipo is localized within different intervals 
depending on the value of gi. These are given by 



-7r/3 < ^0 < 0, 
< V^o < 7r/3, 
vr/3 < V^o < 27r/3, 
L^o = 27r/3, 



^ (1) 

gi<gc- -g2SA, 

(1) / / (2) 

gy < gi<gc' = g2SA 

„(2) ^ „ < „(3) 

gc < gi s gc 



2g2SA + 9/8 s 



.3/2 
'A ' 
3/2 
A ; 



(39) 



gi > g' 



(3) 
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The end points of the intervals in Eq. ()39|). ipQ = kn/S with —l<k<2, represent the 
uniaxially anisotropic structures at the surface 



Sa(0) = sa(3z®z-I)/2, ^o = -vr/3, (40) 

Sa(0) = -sa(3x®x-I)/2, ^0 = 0, (41) 

Sa(0) = SA(3y®y-I)/2, = vr/3, (42) 

Sa(0) = -sa(3z®z-I)/2, ^o = 27r/3. (43) 

From Eqs. and (fl^ it is seen that, for the angles ipo = — vr/3 and ■j/'o = 7r/3, surface 



alignment will be homeotropic and homogeneous (monostable planar), respectively. The 
structure ()41|) coincides with uniaxial ordering in the bulk Q and the surface order param- 
eter tensor ()43p corresponds to planar (random in-plane) alignment. 

According to Eq. ()39|). the case of planar alignment occurs only if the coupling constant 
Qi is positive and the bulk order parameter sa is below its critical value Sc defined by the 
relation 

91 = 2g2Sc + 9/8 s^^'. (44) 

Fig. H shows that the critical order parameter Sc is an increasing function of gi. 

In Fig. Efa) we have plotted the curves representing the components of the order param- 
eter tensor at the surface in relation to the order parameter in the bulk to illustrate that 
destruction of the planar alignment takes place in a second order transition manner. 

However, it should be stressed that our model becomes inapplicable in the immediate 
vicinity of Sc where ^-dependence of the biaxiality parameter Pa critically slows down. Ac- 
tually, as it can be inferred from Fig. El the characteristic length of spatially varying biaxiality 
diverges logarithmically as sa approaches Sc from above. Under these circumstances, the 
assumption that scale of the spatial variations is much shorter than the film thickness is no 
more justified. 

By contrast to the boundary conditions with gi > 0, there are no second order transitions 
provided the coupling constant gi is negative. At sufficiently large values of \gi\, the sur- 
face ordering remains nearly homeotropic. Otherwise, the surface order parameter changes 
smoothly with sa towards the bulk order parameter (jH]). From Eq. (jSHl), for gi = —g2SA, dif- 
ference between the order parameters vanishes. The curves presented in Fig. Efb) illustrate 
this point. 

Leaving aside a detailed discussion of what happen when the coupling constant g2 is 
negative, we just note that in this case the above discussed transition will generally be first 
order leading to jump-like behavior of the order parameter at the surface. 

III. EXPERIMENT 

Now we pass on to describing the experimental procedure employed to obtain the data 
linking the anchoring energy strengths and the dichroic ratio as a measure of the photo- 
induced ordering. To this end we carried out the absorption spectra and the anchoring 
energy measurements for the azo-dye films irradiated at varying exposure time. Thus, we 
used the samples prepared at different irradiation doses to measure the anchoring strengths 
and the dichroic ratio in relation to the dose. The data then can be recalculated to obtain 
the required anchoring energy vs dichroic ratio dependence. 
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FIG. 4: UV-visible absorption spectra of nonirradiated and irradiated azo-dye films. 



A. Sample preparation 



Following the method described in Ref. [21[, the azobenzene sulfuric dye SD-1 was syn- 



thesized from corresponding benzidinedisulfonic acid using azo coupling. The azo-dye com- 
pound SD-1 was mixed with the polymerizable azo-dye SDA-2 in the ratio 40% to 60%. The 
mixture was dissolved in N,N-dimethylformamide (DMF) and a heat initiator V-65 (from 
Wako pure chemical industries, Ltd.) that was added in relation of 1:50 to SDA-2. 

The solution was spin-coated onto glass substrates with indium-tin-oxide (ITO) electrodes 
at 800 rpm for 5 seconds and, subsequently, at 3000 rpm for 30 seconds. The solvent was 
evaporated on a hot plate at 100° C for 10 minutes. 

The surface of the coated film was illuminated with linearly polarized UV light using 
super-high pressure Hg lamp through an interference filter at the wavelength 365 nm. The 
intensity of light irradiated on the film surface at varying time exposure was 2.7 mW/cm^. 
After the photoaligning procedure, the SDA-2 films were polymerized by heating at 150°C 
for 1 hour in vacuum. In order to recover quality of the photoalignment degraded after the 
polymerization, the films were exposed to the UV light for 1 minute regardless of the initial 
time exposure. 

Two glass substrates with the photoaligned films were assembled to form liquid crystal 
cells to measure the azimuthal and polar anchoring energy strengths. The cell thickness 
was 5 /xm and 18 /xm, respectively. Liquid crystal mixtures MLC-6080 (from Merck) in an 
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Irradiation dose (J cm ) Irradiation dose (J cm ) 

FIG. 5: Dependence of (a) absorption coefficients and (b) absorption order parameters on UV 
irradiation dose. Solid lines give interpolation of the experimental data calculated using Akima 
splines. 

isotropic phase were injected into the cell by capillary action. 



B. Absorption spectra 

The UV-visible absorption spectra of the films were measured in the spectral range from 
250 nm to 600 nm for the normally incident probing light which is linearly polarized parallel 
(along the x axis) and perpendicular (along the y axis) to the polarization vector of the 
activating light. 

For nonirradiated films, the curve shown in Fig. |3] as a solid line demonstrates that the 
absorption coefficient does not depend on the polarization state of the testing beam. By 
contrast, as it is illustrated in Fig. EJ the absorption coefficients D\\ and D± differ in the 
irradiated films, thus, revealing the light induced absorption dichroism. This dichroism is 
mainly caused by the photo-induced angular redistribution of the azo-dye molecules. 

By varying the exposure time the films were prepared at different irradiation doses and 
the optical density components and D± at the absorption maximum of azo-dyes (A^, ~ 
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FIG. 6: Azimuthal anchoring energy as a function of (a) UV irradiation dose and (b) dichroic 
ratio ^ which is proportional to the order parameter of azo-dye molecules, 5a (see Eq. Q). Solid 



line represents the result of linear fitting W^f, 



(0) 



w^a^'R with = 0.07225 mJ/m^ and 



(1) 



0.4503 mJ/m^ 



350 nm) were estimated from the measured absorption spectra. The dichroic ratio then 
can be computed from the formula Q. The results for the absorption coefficients and the 
absorption order parameters, which are proportional to the dichroic ratio, are presented in 
Figs. Efa) andEfb), respectively. 



C. Anchoring energy strengths 

The azimuthal ancho ring strength, IV^, was measured in a twisted nematic cell using the 
torque balance method (42 . l43l |. The azo-dye aligning film and a rubbed polyimide layer 
were used as confining substrates. The twist angle was 90 degrees. 

Measurements of the polar anchoring strength. Wo , in anti-parallel aligned cells were 
carried out using the high- voltage technique (22, 13, 0, li^ . The experimental data for the 
azimuthal and polar anchoring strengths are plotted against the irradiation dose in Fig. IHfa) 
and Fig. Efa), respectively. At low irradiation doses, when the exposure energy is below 
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FIG. 7: Polar anchoring energy as a function of (a) UV irradiation dose and (b) dichroic ratio Q 
which is proportional to the order parameter of azo-dye molecules, Sa (see Eq. 0). Solid line 
represents the result of linear fitting —w^^ — w^^ R with w^'^ = 0.38714 mJ/m^ and w^^ = 

3.4116 mJ/m^. 

1 J/cm^, our experimental technique failed to provide accurate estimates for the anchoring 
strengths because of poor quality of NLC alignment in the cells. 

The experimentally measured dependence of the dichroic ratio on the irradiation dose can 
now be combined with the results for W^p and We so as to recalculate the anchoring energy 
strengths as functions of the dichroic ratio, R. The resulting data are shown in Figs. Efb) 
andEfb). 



IV. RESULTS 

As is shown in Fig. ISfa) representing the absorption coefficients, D\\ and D±, measured 
in the film irradiated at various irradiation doses, within the limits of experimental error, 
the average absorption coefficient, Dav, defined by the relation (jH)), remains unchanged at 
irradiation doses higher than 1 J/cm^. So, from the discussion given at the end of Sec. Ill Al 
we conclude that the azo-dye order parameter in the bulk of the film is uniaxial and is of 
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FIG. 8: Azimuthal anchoring energy versus azo-dye order parameter. Solid lines represent the- 
oretical curves computed from the formula (jHUl) by minimizing the energy (jSHJ for two different 
sets of the parameters: (a) gi = 0.47, 52 = 0.4, = 0.067 mJ/m^; (b) gi = —0.9, 52 = 0.4, 
= 0.157 mJ/m^. 



the form given by Eq. 0. In Fig. Efb), it is indicated that the components of the order 
parameter tensor are proportional to the dichroic ratio 0- Clearly, for > -Dy and 
cry > a_i, the azo-dye order parameter Sa and R are both negative. 

The experimental results for the azimuthal and polar anchoring strengths measured in 
NLC cells with photoaligned azo-dye films used as aligning substrates are presented in 
Figs, ini^a) andlZl^a), respectively. The films differ in an amount of photoinduced anisotropy 
which is controlled by varying exposure time and the anchoring strengths are plotted in 
relation to the irradiation dose. 

However, the fundamentally important characteristic describing degree of the photoin- 
duced anisotropy is the azo-dye order parameter. So, in order to compare the experimental 
data and the theory, we need to relate the anchoring strengths and the dichroic ratio. Com- 
bining the anchoring energy data and the curve depicted in Fig. [Sfb) gives the result shown 
in Figs, ini^b) and[7|^b). 

From the relations ()30 |) - (jnnj) the anchoring strengths depend linearly on the dichroic ratio 
provided the order parameters at the surface do not differ from their values in the bulk. The 
results of linear fitting of the experimental data are shown as solid straight lines in Figs. IH^b) 
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FIG. 9: Polar anchoring energy versus azo-dye order parameter. Solid lines represent theoretical 
curves computed from the formula (|H2)) by minimizing the energy (jHEj) for two different sets of the 
parameters: (a) gi = 0.47, 92 = 0.4, = 1.113 mJ/m^, wf^ = 2.711 mJ/m^; (b) gi = -0.9, 
c/2 = 0.4, wP = -0.382 mJ/m2, u,^ ^ 2.247 mJ/m^. 

andlD^b). 

Referring to Fig. Efb), the linear approximation for W^p predicts that the azimuthal 
anchoring strength vanish at certain non-zero value of the dichroic ratio, R ~ —0.16. By 
contrast, from the formulas ()3Up and (jHTj) the anchoring strength is proportional to 
R and, thus, disappear only in the limit of weak photoinduced anisotropy where i? — > 0. 
Assuming that this discrepancy can be attributed to the difference between the bulk and 
surface order parameters of azo-dye, we can apply the phenomenological model described in 
Sec. IIIDI to interpret the experimental data. 

In the angle-amplitude representation (jHS)), the surface order parameters that enter the 
expression (jHUjl are given by 

Si^'^\z=o = -SACOstpo, 5'^y^|2=o = sacos(V'o - 7r/3), (45) 

where —sa = Sa = {^cT/crav)R is the scalar order parameter in the bulk of the azo-dye film 
[see Eq. ®]. According to our model, ipQ is the angle that minimizes the energy (jHHj) . 

So, the computational procedure involves two steps: (a) minimization of the energy ()38p 
to find the angle tpQ] and (b) using the relations to compute the azimuthal anchor- 
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ing energy (IHUj). Following this procedure, we may calculate dependence of the anchoring 
strength on the photoinduced order parameter sa- 

As it is discussed in Sec. Ill Dl the result crucially depend on the boundary conditions that 
are determined by two coupling constants, gi and g2. At gi > 0, the surface favors planar 
(random in-plane) alignment of the azo-dye molecules. In the opposite case of negative 
coupling constant gi, the alignment is homeotropic. 

In figure |Hl we show the theoretical curves calculated for both planar and homeotropic 
boundary conditions. The corresponding numerical results for the polar anchoring strength 
are presented in Fig. El 

The curve plotted in Fig. |Hfa) indicates that, for planar boundary conditions, the az- 
imuthal anchoring strength W^p takes non-zero values and starts growing only if the azo-dye 
order parameter sa exceeds its critical value Sc- Such threshold behavior is a consequence of 
the second order transition discussed in Sec. Ill Dl Contrastingly, as is shown in Fig.|Hl^b), Wif, 
is a smoothly increasing function of sa when the boundary conditions favor the homeotropic 
alignment at the surface. 

At first glance, the curves representing the polar anchoring strength We plotted against 
Sa in Figs, infa)-(b) do not show any noticeable differences. It, however, should be stressed 
that the planar boundary conditions prevent the polar anchoring energy from decaying to 
zero as the order parameter sa decreases. From Fig. Efb) it can be seen that this is no longer 
the case when the coupling constant gi becomes negative. 

These results show that interplay between photoinduced ordering in the bulk of the azo- 
dye films and the preferred alignment of the azo-dye molecules at the surface may have a 
profound effect on the order parameter dependence of the anchoring energies. For the planar 
alignment with gi > 0, our model predicts that the surface ordering may change through the 
second order transition as the photoinduced anisotropy increases. This transition bears close 
resemblance to the second order transitions in nematic liquid crystals previously studied in 
Refs. 0,113. 

Quantitatively, the polar anchoring energy Wg appears to be an order of magnitude higher 
than the azimuthal energy W^. From our estimates, the ratio of the coefficients w^g^ and 
w^, w^q'^ /w^ = —{cf^ + c^^^)/cf^ , is likely to be well above 10. 

The coefficient c^'* is negative and its sign is determined by the dominating contribution 
from the Maier-Saupe term V220 of the spherical harmonics expansion for the azo-dye - NLC 
intermolecular potential Ua-n- From Eq. (jBlOjl the absolute value of c^^ can be significantly 
reduced when the quadrupole term ^224 and the harmonics V222 are predominately positive, 
so that Jg°° zv224,{z) d^; > and Jq°° -2^222 (-z) d-z > 0. 

Under these circumstances, the condition (jHlj), that requires the sum of cf^ and cf^ to be 

(2) (3) 

positive, can be satisfied only if the contribution of the quadrupole term to the sum c)^ +0^ 
is dominating. Thus, we may conclude that the quadrupole term is of vital importance for 
the understanding the reasons behind the significant difference in magnitude between the 
photoinduced parts of the polar and azimuthal anchoring strengths. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we have studied both theoretically and experimentally effects of photo- 
induced ordering in the azo-dye aligning films on the anchoring energy strengths. These 
effects are governed by dependence of the strengths on the azo-dye order parameter. 



18 



Our theoretical approach rehes on the mean field theory j28[ and provides general ex- 
pressions for the Landau-de Gennes surface free energy (fT^ and the anchoring energy ()22p . 
The theoretical results for the azimuthal and polar anchoring energy strengths are obtained 
under certain simplifying assumptions and used to interpret the experimental data relating 
the anchoring strengths and the dichroic ratio. 

We found that linear fitting of the data for the azimuthal anchoring strength, though 
giving good results, predicts, in contradiction to the bare anchoring theory, the effect of 
vanishing anchoring which occurs at certain non-zero value of the dichroic ratio (and, thus, 
the irradiation dose). By using a simple phenomenological model we have shown that this 
effect can be attributed to an interplay between the light induced ordering in the bulk and 
the boundary conditions at the surface of the film which may counteract the action of light. 
Thus, the bulk and surface values of the azo-dye order parameter are generally different. 

For planar boundary conditions that favor the random in-plane alignment of the azo-dye 
molecules, our model predicts threshold behavior of the azimuthal anchoring strength that 
starts growing provided the bulk order parameter, sa, exceeds its critical value. When the 
boundary conditions are homeotropic, this is no longer the case and the azimuthal strength 
smoothly increases with the dichroic ratio. 

The results presented in Figs. |S1 and M demonstrate that the theoretical curves calcu- 
lated for both types of the boundary conditions can fit the experimental data equally well. 
There are, however, differences concerning the polar anchoring strength. By contrast to the 
homeotropic boundary conditions, it never equals zero at the planar conditions. From the 
previously published results li^ it can be concluded that the homeotropic boundary 
conditions is unlikely to occur in the azo-dye films under consideration. 

Our final remark concerns some of the simplifying assumptions taken in our theoretical 
analysis. A more sophisticated theory that goes beyond the Fowler approximation (fTUI) is 
required to take into account surface adsorption phenomena. A self-consistent treatment of 
two order parameter tensors in the interfacial layer also remains a challenge. We hope that 
our results will stimulate further progress in the field. 
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APPENDIX A: IRREDUCIBLE TENSORS, ORDER PARAMETER AND IN- 
VARIANTS 

In this appendix we introduce notations and definitions used throughout the paper. In 
addition, we express the order parameter in terms of irreducible tensors and deduce a number 
of algebraic relations simplifying the derivation of the tensorial form of the intermolecular 
potential given in Appendix iBl 

The irreducible tensors, T^, with the azimuthal number m ranged from —2 to 2 can be 
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defined as linear combinations of tlie following form |50|: 

1 

Tr.= E Cllle.^e^, (Al) 

fl,U= — l 

where C^J^^^ is the Clebsch- Gordon (Wigner) coefficient and e±i = =f(x ± iz)/\^, = z 
are the vectors of spherical basis (x, y and z are the unit vectors directed along the corre- 
sponding coordinate axes). Substituting the values of the Wigner coefficient into Eq. ()A1|) 
gives the expressions for 

To= (3eo®eo-I)/v^, (A2) 

T±i = (eo ® e±i + e±i ® eo)/V2, (A3) 
T±2 = e±i®e±i, (A4) 

so that it is not difficult to verify the validity of the orthogonality relation 

Tr[T„T_„] = {-ir6mn (A5) 

and the algebraic identities 

Z ■ T„ ■ Z = a/273 (5^0, 2 -TrnT-n- ^ = C\m\Smn, (A6) 

where Cq = 2/3, Ci = —1/2 and C2 = 0. 

Under the action of rotation the vectors of spherical basis transform as follows 

1 

e;.^e^(u) = Dl^i^)e,, (A7) 

eo(u) = u, (A8) 
e±i(u) = T(e..(u) ± iey{u))/V2, (A9) 



where D^^(u) = D-'^^{9,(f)) is the Wigner D function j50|, |51|; 9 and are Euler angles of 
the unit vector u; e^(u) = (cos ^ cos 0, cos 6* sin 0, — sin ^), ej^(u) = (— sin 0, cos0, 0). 

The definition ()A1|) implies that transformation properties of the tensors under ro- 
tations are determined by the irreducible representation of the rotation group with j = 2, 
where j is the angular momentum number. So, we have 

1 2 

T„ ^ T^(u) = Cllle,{u) ® e,(u) = ^ ^L(fi)T.. (AlO) 

^i,l'=-l k=-2 

Eq. (jAlOjl can now be combined with the relations ()A8jl and ()A2j) to yield the expression 
for the order parameter tensor Q(u): 

Q(u) = ^372X0(6) = (3u®u-I)/2, (All) 

where the unit vector u is directed along the long molecular axis. 
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The director fi is defined as an eigenvector of tlie orientationally averaged order parameter 
tensor 

(Q(u)), = {DUct>\ 9')) T,(n), (A12) 

k=-2 

wliere 9' and 0' are Euler angles of the vector u related to the basis vectors ej(n). 
Since fi is the director, the averages {|D^^q(0', 6*')) vanish. Other averages 



{DUcI^\9')) = ^S, (A13) 

(^l2o(0',^')> = Pexp(±2^7)/v^ (A14) 

are proportional to the scalar order parameter 5* and the biaxiality parameter P. 

By using the orthogonality conditions ()A5|1 and Eqs. ()A12j) - ()A14jl we recover the relations 
in the traditional form 30j: 

(Q(u)) - = S(n) = S Q(n) + P(m ® m - 1 ® i)/2, (A15) 

^= ^3 (u -11)^-1^/2, (A16) 

P = 3((u-m)'-(u-y\/2, (A17) 



where m = cos7e2;(n) — sin7ey(n) and I = sin7e^(n) + cos7ey(n). 

The order parameter ()A15j) is a traceless symmetric tensor. Therefore, there are two 
non-vanishing independent invariants 

h = Tr[S2(n)] = (3^2 + P^) /2, (A18) 
h = Tr[S=^(n)] = 3S {S^ - P^) /4, (A19) 

which enter the non-elastic part of the well known phenomenological expression for the 
Landau-de Gennes free energy density 

/lg = y (2^ - - ^ ^3 + ^ II (A20) 

where T is the temperature and T* is the supercooling temperature. 

For the scalar order parameters ()A16|1 and ()A17j) combined into a pair {S,P), it is con- 
venient to introduce what might be called the "polar" (or amplitude-angle) representation 

S = scosip, P = s sin?/', (A21) 

where = 2J2/3. Using the representation ()A21|) the free energy density (|X20|l can be 
recast into the form 

kcis, ip) = a{T - T*)s^ - Bs' cos(3^) + = B'C-'U^g, (A22) 
^^LG(//,V') = ^^'-^'cos(3V')+r7^ V = ^s, (A23) 

where t = 32aCB^'^{T — T^) is the dimensionless temperature parameter and = 
T* + B'^/{4aC) is the temperature of the bulk nematic-isotropic transition. The rescaled 
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density ()A23|) is a generalized version of the dimensionless free energy density previously 
used in Refs. ji^, Ii3 • 

Finally, we write down the components of the order parameter tensor ()A16|) in the polar 
representation 

Sij = s [niUj cosip + rriimj cos{ip — 27r/3) + Ulj cos^ip + 27r/3)] (A24) 

and notice that the stationary points of the free energy density ()A22|) where the angle ip is 
a multiple of tt/S represent uniaxially anisotropic states. The latter immediately recovers 
the well known result about uniaxial anisotropy of NLC equilibrium states js^]- 



APPENDIX B: TENSORIAL FORM OF INTERMOLECULAR POTENTIAL 



We begin with the intermolecular potential between two rigid, axially symmetric 
molecules expanded in a series of spherical harmonics as follows ^53>] 



f/(r,Ui,U2) = ^ Mjijy(r) C^XL 

ji,j2,j mi,iri2,iri 

xy,,^,(ui)r,,„ju2)y}:„(f) (bi) 

where r = ri2 = ri — r2, fj and Uj are the position and orientation (equivalently, Euler angles 
of the long molecular axis) coordinates of the interacting molecules, respectively; l^m(u) = 
A/(2j + l)/(47r)i5^*Q(u) is the spherical function [13, 113 • The form of the expansion ()B1|I 
implies that the potential is invariant under translations, — > + Ar, and rotations, 
{rj,Uj} — s> {i?rj,i?Uj}. In addition, we shall assume the head-tail symmetry 

U (r, ui, U2) = U (r, -ui, U2) = U (r, Ui, -U2), (B2) 



so that the outer sum in Eq. (jBlll is restricted to run over even values of ji and j2- 
It is now our task to link the pairwise potential integrated over in-plane coordinates 

l^(z,Ui,U2)= / f/(r,Ui,U2)dxd2/ 
Js 

= J2 [(2ji + 1)(2J2 + l)(2j + l)/(47r)T/' 

X Vnj2A^) Yl ^^o(fii) ^i^„.o(fi2) Dl,{z) (B3) 

m 

and the tensorial representation that was originally suggested by Ronis and Rosenblatt in 
Ref. 

V{z, ui, U2) = V{z, Qi, Q2) = V,Uz) + z ■ (Qi + Q2) ■ z 

+ f3^'\z) Tr(QiQ2) + (3^'\z) z ■ Q1Q2 • z 

+ /5(=^)(^)[z-Qi-z][z-Q2-z], (B4) 

where Qj = Q(uj) is defined by Eq. ()A11|) . In other words, the problem is to express the 
coefficients P^^\z) in terms of the harmonics Vj^j^j{z). To this end we restrict ourselves to 
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the lowest order harmonics of the expansion ()B3|) with ji < 4 and consider the equation 
xi [z ■ To(ui) • z][z ■ To(u2) • z] + X2 Tr(To(ui)To(u2)) 

2 

+ X3 Z ■ To(Ui)To(U2) ■i = Yl «-^mo(fil)^-mo(U2), (B5) 

m=0 

that need to be solved for xi, x^ and X3. The sum on the right hand side of Eq. ()B5|1 
represents sum of the harmonics with ji = = 2. The case of the harmonics with jij2 = 
is much easier to treat as we only have to use the relation 

z-To(u)-z = Do'o(u)- (B6) 

By using Eq. ()A10|) combined with the relations ()A5|) and ()A6|) it is straightforward to 
transform Eq. ()B5|) into a system of linear equations. The solution of the system is given by 

xi = (3ao + 4ai + a2)/2, 
X2 = -2{ai + 0:2), 

X3 = a2. (B7) 

Given the values of the Wigner coefficients, we can now use the relations ()B6|) and ()B7|) 
to derive the final result in the following form: 

VUz)=Vooo{z)/i^7if\ (B8) 
[(4vr)3/V5] = V202{z) = vo22{z), (B9) 

+ 10/VUV222{Z) + 3/VUv224{z)}, (BIO) 



[(47r)3/V5]/5(2)(z) = -20/V70{v222{z)+V22a{z)}, (Bll) 



[(47r)3/V5] (3^'\z) = 35/v^t;224(^). (B12) 

The formulas ()B8|) - ()B12|) relate the parameters of the representation ()B4|) to the coeffi- 
cient functions in the spherical harmonics expansion ()B3jl . The terms f202 and ^022 describe 
the coupling between orientation of the molecules and the intermolecular vector, whereas 
f22o and f224 are known as the Maier-Saupe and the quadrupole terms, respectively. For 
the interaction between NLC molecules, the functions P^^\z) and f3^'^\z) define the elastic 
coefficients of NLC that must be positive. This stability condition implies that P^^^ and 



are both predominately non-positive |27l l28.. .31. .SSi]. 
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